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This article presents a derivation of the transport mean free path for magneto-transverse light 
diffusion, £*^, in an arbitrary random mixture of Faraday-active and non- Faraday active Mie scat- 
terers. This derivation is based on the standard radiative transfer equation. The expression of the 
transport mean free path obtained previously from the Bethe-Salpether equation, for the case where 
only Faraday-active scatterers are present, is recovered. This simpler formulation can include the 
case of homogeneous mixtures of Faraday-active and non-Faraday active scatterers. 



I. INTRODUCTION 



Magneto-transverse light diffusion - also known as the "Photonic Hall Effect" (PHE) - was theoretically predicted 
five years ago by Van Tiggelen ||l|], and was experimentally confirmed one year later by Rikken This effect 
is analogous to the electronic Hall effect, well known in semiconductors physics. The evident driver behind the 
electronic Hall effect is the Lorentz force acting on charged particles. The PHE finds its origin in the Faraday effect, 
present inside the dielectric scatterers, which slightly changes their scattering amplitude. This is the reason for the 
similarities between the PHE and the so-called Beenakker-Senftleben effect, which concerns the transport coefficients 
of dilute paramagnetic gazes 

Motivated by the experimental observation of the PHE, theoretical work was first started on the single scattering 
of spherical magneto-optical particles. The scattering matrix and the scattering cross-section were calculated exactly 
for a single Faraday-active dielectric sphere using perturbation theory Q. From this solution, the Stokes parameters, 
which completely describe the intensity and polarization of the scattered light, were derived Perturbational 
methods for the scattering by weakly anisotropic particles were first used by Kuzmin et al. for the single scattering 
case 

The PHE in multiple scattering is controlled by a length, which has been defined as the mean free path for 
magneto-transverse light diffusion, i*^. A simple expression for this length was obtained and successfully compared 
to experiments using a formulation based on the ladder approximation of the Bethe-Salpether equation [Q. The 
experiments investigated the dependence of the PHE on the volume fraction of the scatterers, first on the real part 
of the dielectric constant of the scatterers ||^, and more recently on their imaginary part However, some parts of 
the derivation of Ref. are rather technical and fail to give a satisfying explanation of the origin of the dependence 
of t*^ on the anisotropy of the scattering. This article presents a simpler derivation of t*^, based on the radiative 
transfer equation, which should clarify this point. The radiative transfer equation is a Boltzmann type equation 
which describes the transport of light in multiple light scattering From the radiative transfer equation, the 

diffusion equation is derived in an infinite medium when the scattering is not highly anisotropic. In other types of 
anisotropic materials, such as single-domain nematic liquid crystals, similar methods were developed, using either the 
Bethe-Salpether equation [0 or the radiative transfer equation Q . 

II. SINGLE SCATTERING 

The single scattering of light by one dielectric sphere made of a Faraday-active material embedded in an isotropic 
medium with no magneto-optical properties is first considered. In a magnetic field, the dielectric constant of the 
sphere, eb, is a tensor of rank two. It depends on the distance to the center of the sphere of radius, i?, via the 
Heaviside function, 6(|r| — i?), which equals 1 inside the sphere and outside 

eB(B, r)jj = (eo - l)5ij + lepeijkBk 6(|r| - R), (1) 

where eo is the value of the normal isotropic dielectric constant of the sphere, ep = 2eo^^^Vo-B/w is the coupling 
constant of the Faraday effect, 6 and e are respectively the Kronecker and Levi-Civitta tensors. The Verdet constant 
of the Faraday effect is denoted by Vb, -B is the amplitude of the magnetic field and uj is the frequency. The intensity 
of the scattered light, in single scattering, is characterized by the phase function which is also the essential ingredient 
for the transport theory of light. The phase function, F(k, k', B), is proportional to the differential scattering cross- 
section averaged with respect to incoming and outgoing polarizations. It depends on the direction of the incoming 
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plane wave, k, on the direction of the scattered wave, k', and on the direction of the magnetic field, B. The hat above 
the vectors denotes normalized vectors. For spherical magneto-optical particles and to linear order in the applied 
magnetic field, this phase function can be written as Q 

F(k,k',B) =^^o(k,k')+det(k,k',B)Fi(k,k'), (2) 

where det(A, B, C) = A • (B x C) denotes the scalar determinant. Due to the rotational symmetry of the scatterer, 
the functions, Fq and Fi, only depend on the scattering angle, 0, which is the angle between k and k'. The phase 
function Fi only depends on the difference in the azimuthal angles associated with k and k', because of the axial 
symmetry of the scatterer around the direction of the magnetic field. These symmetry properties simplify considerably 
the use of the radiative transfer equation [|l3| . The amplitude of Fi was found to be proportional to the dimensionless 
parameter, ep 

The albedo is introduced as the ratio of the scattering cross-section over the extinction cross-section jlj] 

Qscatt /o\ 

a=-^ ■ (3) 

It is related to the phase function defined in Eq. (|^), 



a 



= / dkF(k,k',B) = / dkFo(k,k'), (4) 



where the last equality follows from Eq. (0). Thus in the absence of absorption the phase function is normalized 
{a=l). 



III. TRANSPORT IN MAGNETO-OPTICAL DIFFUSE MEDIA 



In an isotropic medium much larger than the transport mean free path, £* , Fick's law relates the diffusion current, 
J, to the energy-density gradient, V/, by J = —DqSIF The conventional diffusion constant for radiative transfer 
equation is denoted Dq and is usually related to the transport mean free path, i* , and the transport velocity, wg, 
by Dq — ^ve(*- In the presence of a magnetic field, the diffusion constant is a second-rank tensor. The part of 
this tensor linear in the magnetic field is responsible for the PHE, whereas the part quadratic in the magnetic field 
generates a photonic magneto-resistance, which could also be observed experimentally By Onsager's relation, 

Dij(B) — Z?ji(— B), the part linear in the external magnetic field must be an antisymmetric tensor, and Fick's law 
becomes, 

J = -D(B) • V/ = -DoV/- L'^lB X V/. (5) 

The term containing D± describes the magneto-transverse diffusion current responsible for the PHE. In analogy to 
the definition of £* , the transport mean free path, £*j_, for magneto-transverse diffusion is defined as D± — ^ve£*i- 
For the electronic Hall effect, is proportional to the Hall conductivity, axy, the sign of which is determined by the 
charge of the current carriers. Similarly, a positive means that the PHE has the same sign as the Verdet constant 
of the scatterers, and a negative £*j_ means that the PHE has an opposite sign, which is also possible depending on 
the scattering. This point has been carefully checked experimentally ||^. 

The phase function introduced above can be used to describe not only the single scattering of light by independently 
scattering particles but also light scattering by a large collection of such particles in the regime of multiple light 
scattering. It is assumed that only magneto-optical Mie scatterers are present, and included in a matrix with no 
magneto-optical properties. A transport equation for a medium comprising randomly distributed spherical particles 
embedded in an isotropic medium can be written as follows [ p^ : 

— dtI(r,k,t)+ek-VJ{r,k,t)+I(r,k,t)^ / dk' i^(k, k', B)I(r, k', i), (6) 

VE J 

where £ denotes the elastic mean free path, which is the average distance between two subsequent scattering events. 
The specific intensity Z(r, k, t) is defined as the density of radiation at position r, time t, in the direction k. The 
specific inten sity has the dimension of a monochromatic radiance: energy per unit solid angle, time, wavelength and 
surface area p7|]. When supplemented by appropriate boundary conditions, the radiative transfer equation can be 



2 



solved for a particular problem. It is useful to introduce the average specific intensity, /(r, t), and the current, J(r, t), 
which are defined as 

I{r,t) = J J(r,k,t)dk, J{r,t) =ve j kdkl(r,k,i). (7) 

To simplify the notations, the dependence of the specific intensity or of related quantities upon the magnetic field, B, 
is not explicitly indicated. The integration of Eq. (|^) with respect to k leads to the continuity equation 

at/(r,<) + V.J(r,t) = -il^^/(r,i), (8) 

where a denotes the albedo defined in Eq. (^). 

If Eq. (||) is multiplied by k and integrated over the direction k, the right-hand side of Eq. (||) will contain two 
integrals which depend on and Fi. These integrals are the two dimensionless quantities: 



< cos6' >=y dki^o(k,k')k-k', (9) 

and 

Ai =y"dkFi(k,k')det(k,k',B)2. (10) 

By choosing a system of coordinates linked with k', it can be proved quite generally that 

dki^(k, k', B) k =< cos 61 > k' - AiB X k'. (11) 



After the integration over k' the following equation is obtained 

VE 



dt3 + 3 (1- < cos 61 >) + AiB X J = -IveV ■ j dkX(r, k, t)kk. (12) 



If it is assumed that the angular distribution of the specific intensity is almost isotropic, the current is much smaller 
than the average specific intensity. In that case the following approximation is valid: ||l^ 

J(r, k, t) « /(r, <) + — J(r, t) ■ k + .... (13) 

VE 



Substituting this relation into Eq. (|12[), and neglecting dt3 for processes which are slow with respect to the charac- 
teristic time between two scatterings, I/ve, gives Fick's law Q: 

J(r,i) = -D.V/(r,i). (14) 

The diffusion tensor is given by, 

D(B),y = \vEi [(1- < cobO >) 5,^ - A^e,,kBk] ' • (15) 



It is important to note that the magnetic correction to the scattering cross-section plays a role similar to the role of 
the asymmetry factor, < cos 9 >, present when no magnetic field is applied. The dependence of < cos 9 > can be 
obtained from standard Mie theory [Q. Both anisotropy factors, < cos 9 > and Ai, are non-zero provided that the 
symmetry between forward and backward scattering is broken, which is the case for a finite size scatterer. This is the 
reason why the PHE of a Rayleigh scatterer vanishes whereas the PHE of a Mie scatterer does not Q . 
In the case where Ai <C (1— < cos6' >), an expansion of the bracket gives 

1^1 £ - 

Dij(B) = -VE- ^ — Sij + -veAi- -nSi-jkBk- (16) 

3 l-<cos6'> 3 (l-<cos6i>)^ 

This equation identifies C = ^/(l— < cos9 >) with the transport mean free path, and l*^ — Ail* /{I— < cos9 >) as 
the transport mean free path for magneto-transverse diffusion. The final result of Ref. m is thus recovered with little 
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efforts. The order of magnitude of the ratio P^ji* is 10~^ for an apphed magnetic field of IT and for samples made 
of rare-earth materials 0. The presence in Eq. ( |l6| ) of the coefficient, A\, which is directly connected with the PHE 
of a single scatterer Q , states that the PHE in multiple scattering is directly proportional to the normalized PHE of 
a single Mie sphere, is of the same sign and amplified by the factor, 1/(1— < cos0 >)^. The expansion to first order 
in the magnetic field involved in Eq. (|l6|), clarifies the origin of the factor, 1/(1— < cos0 >)^, which makes l*^ even 
more dependent on the asymmetry factor, < cos^ > than the transport mean free path, . 

Substituting Eq. (|lj) into the continuity equation (^) yields the diffusion equation for the density of radiation 



at/(r,i) = Z?oV2/(r,i) 



(1 - a)vE 
I 



Iir,t)^Do 



V2/(r,t) 



i2 



/(r,t) 



(17) 



As expected, the part of the diffusion tensor, which is linear in the magnetic field and therefore antisymmetric does 
not enter the diffusion equation which depends only on the symmetric part, Dq. In addition to the transport mean 
free path, £*, Eq. ( p^ ) defines the characteristic length for the absorption in multiple light scattering. La 



La = 



I £* -e 

3(1 - a) 



i* ■ L. 



bs 



(18) 



as function of the characteristic length for the absorption of the coherent beam iabs ■ Therefore there is no dependence 



of La on the magnetic field to linear order. As discussed in Ref. 
the absorption cross-section, which is proportional to 1 — a. Eq. 
the medium is dissipative. 



I , the diffusion coefficient Dq does not depend on 
1^ ) shows that this property also holds for i*^ when 



IV. TRANSPORT IN MIXTURES OF FARADAY-ACTIVE AND NON FARADAY-ACTIVE SPHERES 



The main result of Eq. (|16|) can be extended to the case of an homogeneous mixture of Faraday-active spheres with 
volume fraction, uq, and non Faraday-active spheres with volume fraction, ni. The phase function of Eq. (|^) is now 
modified for the case of the mixtures as follows 



F(k,k',B) 



no 



m + ni 
Fo(k,k') 



i^o(k,k') 



"1 



The equation (|T^ can be generalized to 



no + ni 



f^o(k,k')+det(k,k',B)Fi(k,k') 



ni 



m + ni 



det(k,k',B)Fi(k,k'). 



(19) 



I?(B)y = IvEi 



(1- < cos 61 >)Si 



ni 



m + ni 



AiEijkBk 



(20) 



In the case when the last term of the bracket is small with respect to 1— < cos6' >, the following expression for the 
ratio of the transport mean free path for magneto-transverse diffusion over the transport mean free path is obtained 



ni 



no + ni 1— < cos 6 > 



(21) 



V. CONCLUSION 



In conclusion, the present transport theory for magneto-optical spheres of arbitrary size confirms the results of 
Ref. [Q concerning the magneto-transverse light diffusion obtained from the Bethe-Salpether equation. This new 
approach clarifies the origin of the dependence of the transport mean free path for magneto-transverse diffusion on 
the anisotropy of the scattering. It also shows that the magnetic anisotropy in the scattering cross-section plays a 
similar role as (cos 9) in the transport mean free path. This derivation can be generalized to the case of mixtures of 
Faraday-active and non Faraday-active scatterers. 
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